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Abstract 

We present a general Markovian framework for order book modeling. Through our ap¬ 
proach, we aim at providing a tool enabling to get a better understanding of the price 
formation process and of the link between microscopic and macroscopic features of fi¬ 
nancial assets. To do so, we propose a new method of order book representation, and 
decompose the problem of order book modeling into two sub-problems: dynamics of a 
continuous-time double auction system with a fixed reference price; interactions between 
the double auction system and the reference price movements. State dependency is in¬ 
cluded in our framework by allowing the order flow intensities to depend on the order book 
state. Furthermore, contrary to most existing models, the impact of the order book up¬ 
dates on the reference price dynamics is not assumed to be instantaneous. We first prove 
that under general assumptions, our system is ergodic. Then we deduce the convergence 
towards a Brownian motion of the rescaled price process. 


1 Introduction 


Nowadays, most financial exchanges nse a limit order book (LOB) mechanism. In these order- 
driven markets, market participants send their buy and sell orders to a continuous-time double 
auction system, where orders are matched according to their price and time priority. Under¬ 
standing the LOB dynamics is one of the fundamental issues in the held of market microstruc¬ 
ture and leads to many interesting applications in optimal execution, design of electronic trading 
algorithms, minimization of market impact costs, short-term predictions and regulation. In the 
recent years, many works have been devoted to the description of order book dynamics. Order 
book models can be essentially divided into two types: economic models, where one focuses on 


the behaviors of individual agents and their optimal decisions, see for example Parlour (1998), 


Foucault 

(|1999|) and 

Ro§u 

(2009 


and Ro§u (2009); statistical models, where the order hows are seen as random 


processes, see 

Smith, Farmer, Gillemot, and Krishnamurthy 

(2003) 

Cont, Stoikov, and Talreja 

O 

T—1 

o 

CM 

Abergel and Jedidi 

(2011 

), 

Cont and De Larrard 1 

2013 

), 

Lakner, Reed, and Stoikov 


(20131), Lachapelle, Lasry, Lehalle, and Lions (2014), Bayer, Horst, and Qiu (2015) and Abergel 


and Jedidi (2015). With the notable exception of Abergel and Jedidi (2015), where the au¬ 


thors consider the case of Hawkes-type dynamics, these models usually assume Poisson hows 
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for the order arrival processes. Such assumption is mainly made for technical reasons, since 
it is well-known that it is not consistent with market data. In [Huang, Lehalle, and Rosen¬ 


baum (2015), the authors propose to replace the Poisson assumption by a state-dependent 


approach where the intensities of the flows depend on the state of the LOB. This model, called 
“Queue-reactive” model, provides new insights for the order book dynamics, such as market 
participants behaviors conditional on different states of the order book, the LOB’s asymptotic 
form and the bid-ask spread distribution. It is also a very relevant tool for practitioners in the 
perspective of transaction cost analysis of complex trading algorithms. In this paper, we aim at 
extending the Queue-reactive model to a more general framework, in which most of the existing 
statistical models can be included (up to minor modifications). Our goal is to give some theoret¬ 
ical results on the system’s ergodicity as well as the asymptotic scaling limit of the price process. 

In the LOB, price levels are discretized by a minimum price change unity called the tick value 
(denoted by a). Market participants can place their buy/sell orders at any level which is a 
multiple of the tick value and these orders will either stay in the LOB (a buy order with price 
lower than the current best ask price, or a sell order with price higher than the current best bid 
price, this type of orders being called “limit order”), or be matched with the existing orders 
in the LOB (this type of orders being called “market order”). The LOB, as its name suggests, 
is composed of all unmatched limit orders and can be seen as a (rough) approximation of the 
current microstructural supply and demand on the different price levels. 


Current statistical models differ in their way of representing the LOB. In |Cont, Stoikov, and 
Talreja (2010), the price grid is supposed to be hnite {uminOi, ...,nmaxCt), and the LOB is repre¬ 
sented by a Umax —n-min + l dimensional vector that records the buying/selling quantities at each 
of these price levels. In such representation, the different limits are indexed by their absolute 
price level. In practice, to cover the intra-day price range [pmmPmax], the dimension of the state 
space have to be at least -|- ij^ which is typically a very large number. Another way of 

representing the order book state is to use the relative indexing method. Following ideas from 


the Zero-intelligence model of Smith, Farmer, Gillemot, and Krishnamurthy (2003), Abergel 


and Jedidi (2011) propose to use the best bid and the best ask prices as two reference prices 


to index the limits. In that case, the LOB is made of the following elements: the two reference 
prices Pbestud and Pbestask, and the limits around them which are two K dimensional vectors a 
(for the ask side) and b (for the bid side). The vector a = [ui,..., ax] records the limit sizes at 
the price levels [pbestbid + ct, ■■■,Pbestbid + Ka], while the vector b = [6i,..., bx] records the limit 
sizes at the price levels \pbestask — ■■■■,Pbestask — In practice, observing a market depth of 

hve ticks is usually considered enough for most trading purposes. Consequently, for a typical 
stock with spread size of order five ticks, the value of K should be generally of order 10 so that 
essential information from the LOB is captured. Thus, the use of these two reference prices 
reduces signihcantly the dimensionality of the state space. Note that in this representation, the 
index of a limit at a given price level is no longer constant. Therefore, appropriate boundary 
conditions must be defined to deal with price changes. 

In this paper, we propose an original representation of the order book, using only one reference 


^Of course if we choose rimin = 1 and Umax = oo, we have a complete description of all the available 
buying/selling offers in the whole price grid aN. 

^Note that in this representation, there can be overlaps between the vector a and b when the bid-ask spread 
is larger than one tick. 
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price which is not necessarily directly observable from the order book state. We view this ref¬ 
erence price as sort of market consensus about the underlying “efficient” price used by market 
participants when making their trading decisions. We keep K limits on each side of the reference 
price and the LOB is fully described by a 2K + 1 dimensional vector, which is then modeled 
by a continuous-time Markov jump process. The use of this unique reference price gives us a 
lot of flexibility when modeling the order book. Since the reference price is no longer directly 
determined by the order book state, we can differentiate two types of jumps in the Markov 
process: pure order book state jumps (jumps for which the reference price stays constant) and 
common jumps (jumps for which a reference price change is involved). For large tick asset^ 
such decomposition is proved to be very relevant when studying the conditional dependences 
between the dynamics of the LOB and its state, see Huang, Lehalle, and Rosenbaum|| ( |2015[ ). 
Moreover, in this framework, we are able to easily incorporate exogenous price movements into 
the order book dynamics. This can be simply done adding a reference price jump component 
which is independent of the order book state. 


At the high frequency scale, the LOB state is one of the two public information that are accessi¬ 
ble to traders and their automates (the other being the history of the order flows). Thus it plays 
a very important role in their trading decisions. In our framework , the LOB is assumed to be 
a continuous-time Markov jump process, and the influence of the LOB state on the incoming 
flows is modeled through a state-dependent infinitesimal generator matrix for the jump process. 
Indeed, in practice, traders essentially rely on information deduced from the current LOB state 
when deciding to send an order at a specific price level. Various simplifying assumptions on 
the information set used by traders can be considered in our framework in order to facilitate 
the empirical studies. The index of a limit, for example, is probably one of the most important 
elements in their decision process, as it gives the distance between the target price and the 
reference price. Influence of other variables, such as the target limit’s size, its relative distance 


to the current best offer queues and the size of its opposite queue is studied in Huang, Lehalle, 


and Rosenbaum (2015), and are shown to also have non-negligible effects on the dynamics of 


the order flows. 


Under appropriate assumptions, the Queue-reactive model can be easily estimated using em¬ 
pirical data. It provides many interesting new insights on the origin of some micro-structural 
properties, such as the stylized empirical distribution of the LOB state. It has been shown in 


Bouchaud, Mezard, and Potters (2002) that there exist some regularities in the order book’s 
empirical form, that is the average value of the LOB state (a 2K dimensional vector in our 
model). From a theoretical point of view, these regularities are closely linked with the notion 
of ergodicity of the LOB system (the exact definition of ergodicity will be given in Section]^: 
Ergodicity ensures the convergence of the LOB state distribution towards an invariant prob¬ 
ability measure. Thus the stylized form observed on market data might be explained by a 
law of large number type phenomenon for this invariant distribution. This hypothesis is sup¬ 


ported by empirical studies in Huang, Lehalle, and Rosenbaum (2015), in which the authors 


compare the theoretical asymptotic distributions in our model with empirical estimations, and 
show that they are very close. In Huang, Lehalle, and Rosenbaum (2015), some assumptions 


are made to ensure the LOB system’s ergodicity. In this paper, we want to generalize them 
in an extended framework where the volume of the orders is no longer constant and the in¬ 
fluence of the order book state on the dynamics of the reference price may not be instantaneous. 


large tick asset is defined as an asset whose bid-ask spread is almost always equal to one tick. 
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Another important element in order book modeling is the asymptotic behavior of the price. 
Such analysis is very relevant as it provides useful insights on the price formation process, and 
links the dynamics at the microscopic level with macroscopic features of the asset, such as its 
volatility. We prove that in our framework, the rescaled reference price process converges to 
a Brownian motion. An expression for the macroscopic volatility in terms of the flow rates is 
derived using a functional central limit theorem together with the strong mixing property of 
the price increments, in the spirit of Abergel and Jedidi (2011). 


The paper is organized as follows. In Section 2, we set up the general Markovian framework. 
Section 3 discusses the ergodicity properties of the model. The diffusive limit of the rescaled 
price process is stated in Section 4. Finally, in Section 5, we give some specihc models which can 
be seen as particular cases of our framework. The technical proofs are gathered in an appendix. 


2 A general Markovian framework 

2.1 Representation of the order book 

In our framework, the order book is made of two elements: the center position: a certain ref¬ 
erence price (note that here the center position is not necessarily the mid price) and the shape 
of the book (the queue sizes around the reference price), see Figure for an example. The 
center position, denoted by Pref £ {na + 0.5a, n G Z}|^ can be seen as the current consen¬ 
sus price level and is used to index the limits. We write Q±i for the queue at the price level 
Pref (* — 0.5)a and denote its size by q±i and its price by pi. Then the total order book’s 
shape at time t is an infinite vector q(t) = [...,q^k(t), ■■■<l-i(t),qi(t), ...,qk(t),...], with g* G Z, 
and \qi\ is the number of orders at the limit Qi (g* < 0 if these orders are bid orders and g* > 0 
if they are ask orders. Note that in such a representation, one can have g* > 0 or g* < 0 for 
all i). The LOB information at time t is therefore fully represented by X{t) := {q{t),prefit)), 
t > 0. To restrict the value of X{t), we consider only K limits on each side, and thus have now 
q{t) = [g-ic(t), ...g-i(t), gi(t),..., qxit)] and X{t) G x a(0.5 -h Z). 

The state space of q{t) is actually smaller than More specifically, let us dehne 

ibestbidiq) = max{-K - l,sup{i|gi < 0}) 
hestaskiq) = min(A'-h l,inf{*|gi > 0}). 

The state space O of q(t) is defined as the set of all g G such that: for all i G {—K, 
qi <0 if i < ibestbid{q) and g* > 0 if i > ibestask{q)- 


2.2 Dynamics of the order book 

In our general framework, we model the LOB vector X{t) as a continuous-time Markov jump 
process, whose infinitesimal generator matrix Q will be given in Equation (|^. We differentiate 
two types of jumps in the order book dynamics: pure order book state jumps (for which p^e/ 
stays constant) and common jumps (for which a change in the value of pref is involved). 

^For the generality of the framework, we allow for negative prices, however, in practice, the model should 
of course be used over a reasonably small time interval so that prices remain positive. 
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Figure 1: Limit order book 




Pure order book state jump 

There are thre^ types of orders that interact directly with the order book and trigger pure 
order book state jumps: 

• Limit orders: insertion of a new order in the order book (a buy order at a lower price 
than the best ask price, or a sell order at a higher price than the best bid price). 

• Cancellation orders: cancellation of an already existing order in the order book. 

• Market orders: consumption of available liquidity (a buy or sell order at the best available 
price). 


In the seminal work of Smith, Farmer, Gillemot, and Krishnamurthy (2003), the arrival times of 
the above three types of orders at different price levels are assumed to be mutually independent 
and exponentially distributed. Furthermore, each order has unit size. In our approach, the size 
of the jumps, which represents the amount of volume inserted to/removed from the LOB for a 
given event, is random. Moreover, the arrival rate of a given jump is assumed to be function of 
the index of the target price, the current LOB state vector q{t), the direction of the jump and its 
size. That is, for any q,q' E Q {q q'), p G a(0.5 + Z), n G N"*" and any e* = {a-x, 

{ttj = 0 for j 7^ i and Oj = 1), we have in cases where q + ncj G O and q — nCi G O: 


Q{Q,p),{q-nei,p) 


0, otherwise. 


where the /j and Qi are 2K functions: O x N"*" —?• M"*". 


( 1 ) 


®Four if we also consider the modification orders. 


We view modification orders as a combination of a 


cancellation and an insertion order that arrive in a very short time interval. 
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Note that in Q, fi(q, n) and gi{q^ n) have different meanings for different i and q. For example, 
when i > ibestask, represents the arrival rate of sell limit orders of size n, and gi{q,n) 

the sum of the rate of cancellations of size n and the arrival rate of market buy orders of size 
n. When i < ibestbid, the role of fi{q,n) and gi{q,n) are switched. Note also that q ± nCj is not 
always in the state space 17 even when g G 17. Thus some values of the functions /* and g* are 
not needed in Equation Q and assumed to be equal to zero. Furthermore, so that there is no 
absorbing state, we assume 


+ 3i{Q, n)) > 0. 

i n 


Common jumps 

The reference price Pref can be viewed as a consensus value on the “efficient” price and takes 
discretized values in a(0.5 + Z). In practice, this reference price is built based on two sets of 
information: the current state of the LOB and the historical order flows. We have in mind 
that Pref moves in a Markovian manner, so its dynamics depends on the present information 
only, that is the current state of the LOB. In our framework, we restrict the price jump size to 
one tick at each time. We use two functions m, d : 17 —)■ M+ to describe respectively the rate of 
positive and negative jumps: 


2(9>p)>(g'>p+o) 

= u{q) 

g'eo 


2(9>p)>(g'>P-a) 

= d{q) 

(?'eo 


2(g,p),(g',p±na) 

= 0 , for n > 2 


g'eo 


( 2 ) 


To understand Equation (|^, let us first consider the following simple example where the LOB 
state information is summarized by the difference between the current value of Pref and the 
mid price p^id- 


Example 1. 


u{q) = 
d{q) = 

with Oq > 0 and 9i a positive constant representing the intensity of the adjustment of Pref 
towards pmid- 

In the above example, we assume that the reference price jump rate depends on the deviation 
of the current value of Pref from Pmid- Indeed, Pmid is often considered an approximation of the 
LOB center implied by its current state. We may also include other LOB information such as 
the available quantities at and when defining u{q) and d{q). Such additional 

variables increase the complexity of our model but make it more realistic. 


i^^bestbid T ibestask) /: 

9o + 6*1 max(0, a(imid - 0.5)) 
do + 9i max(0, -a{imid + 0.5)), 
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Remark 2.1. In the already mentioned Queue-reactive model, changes of pref are triggered 
by the order book events that modify Pmid, while in Equation they are driven by the order 
book state. Although the Queue-reactive model cannot exactly be seen as a particular case of 
the framework presented in this paper (see Section^, most of the theoretical results shown in 
Sections 3 and 4 still hold (with some minor modifications in some assumptions). In particular, 
one can prove the diffusive limit of the reference price in the Queue-reactive model applying 
almost the same method as that used here. 

As soon as Pref changes, the value of qi switches immediately to the value of one of its neighbors 
(right if Pref increases, left if it decreases). As we keep only K limits on each side, two bound¬ 
ary distributions ti-k and tik are introduced for generating the new queue sizes at Q-k (when 
Pref decreases) and Qk (when p^e/ increases). To possibly incorporate external information, 
we moreover assume that with probability the LOB state vector q{f) is redrawn from 

, if Pref decreases) when pref changes. As shown in 
), models where price dynamics are purely endogenous, 
driven by order flows only, are usually not able to reproduce some of the important macroscopic 
features of prices, such as the volatility. Thus the parameter 0 ^®*”'** can be understood as the 
percentage of price changes due to exogenous information, in which case market participants 
readjust very quickly their order flows around the new reference price, as if a new state of the 
LOB was drawn from its (invariant) distribution (ergodicity conditions are discussed in the 
next section). 

For g e O, we write g+ = [g.^+i,..., g_i, gi,..., gi^], g" = [g_ic, g-i, gi, gx-i], [q^ fi] = 
[g_i^+i,...,g_i,gi,...,gi^,/] and [l,q~] = [/, g_K, •••, g-i, gi, •••, Under the above assump¬ 

tions, we have for Z G Z and g, g', g" G O such that q’^ 7 ^ g"*" and q''~ 7 ^ g“: 


some distribution (tt®"''® if p^e/ increases 


Huang, Lehalle, and Rosenbaum (2015 


2('?>P)>([9+>Fp+“) 
Q(q,p),W ,p+a) 

Q-{q,p),{[l,q-],p-oi) 

Qiq,p),(<i",p-a) 


(1 - 0®'“®®'*)n(g)7r;^(Z) + 0®’“®®®‘M(g)7r'®®"([g+, /]) 
0®'“®®®*M(g)7r®®®®=(g') 

(1 - 0’'“™*)d(g)7r_ic(/) + 0’'“®®'‘d(g)7r'^'^"([Z, g"]) 
0’'“®®'*d(g)7r‘^""(g"). 


(3) 


The infinitesimal generator matrix of X{t) 

Equations ([^, ([^ and (|^ give a complete description of the inhnitesimal generator matrix Q 
of the process X{t), which is summarized in the following assumption. 

Assumption 2.1. Let q,q',q”,q G Ll, p,p G a(0.5 -|- Z), n G N"*", I G Z be such that g'+ 7 ^ q~^ 
and q"~ 7 ^ q~. The process X{t) is an irreducible Markov jump process with aperiodic embedded 
chain whose infinitesimal generator matrix Q is of the following form (with 2K functions fi,gi : 
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X N"*" —)■ M"*" and two functions u,d : ^ ; 


2(9,p)>(<?+«ei,p) 
Q(q,p),{q -nei,p) 
2(q,p),([g+,i],p+o) 
2 (g,p), ('?',?+«) 

2 (9,P) >(<?",P-«) 
2{g.p),(g,p) 

Q{q,p),{q,p) 


(1 - /]) 

(1 - g-]) 





(g,:p)er2xa(0.5+Z),(g,p)/(q,p) 


0, otherwise. 


(4) 


Note that under Assumption (2.1), the dynamics of the process X{t) is invariant under trans¬ 
lations of the LOB center position: its inhnitesimal matrix generator Q satishes: 


Q(g'^,p^),(q'^,p^+d) ~ 2(ghp^).(9^.p^+/3)’ 

for any q^,q‘^ G O, p^,p'^ G q:(0.5 + Z) and ft G aZ. One can also remark that in our framework, 
the order book state process q{t) alone is also a continuous-time Markov jump process, whose 
ergodicity is discussed in the next section. 


2.3 Comparison with existing models 

The hrst major difference between our approach and the existing Markovian models in the 
literature is the introduction of state dependency in the order book dynamics. Most of the 
current order book models follow the “Zero-intelligence” framework, using Poisson flows for 
the processes of order arrivals. The Poisson assumption is clearly unrealistic, see for exam¬ 
ple the empirical results in Huang, Lehalle, and Rosenbaum (2015). In our framework, we 


propose to incorporate the strategic behaviors of market participants via a mean-held game 
approach, assuming their decisions depend on an underlying “efficient” price pref and on the 
LOB state vector q(t). Note also that Equation ([^ allows us to have jumps of random size in 
the order book’s shape, while a constant jump size is often assumed in the other existing models. 

We also introduce a new method of LOB representation using one unique reference price Pref- 
Most models use the best bid and best ask prices as two reference prices for indexing the buy 
and sell limits. In such models, these reference prices are directly determined by the order 
book state. In particular, changes in the order book state are immediately carried on the 
values of these prices. In our framework, p^e/ is not necessarily deduced from the order book 
state. Therefore, we can assume that changes in the order book state affect the value of pref 
with some delay rate (the functions u and d introduced in Equation ([^). Thanks to this 
original representation, we can naturally decompose the order book dynamics into two parts: 
a continuous-time multidimensional queuing system (Equation ([^) and the dynamics of its 
center, that is the ref erence price (Equations ([2| and ([^ ). Compared with the Queue-reactive 
model introduced in Huang, Lehalle, and Rosenbaum ( 2015[ ), pref is no longer constrained 
within the bid-ask spread. This desirable feature gives us the possibility of separating the 
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exogenous and endogenous parts in the price dynamics by choosing appropriate price jump rate 
functions u and d. For example, with the functions u and d of Example [T| one can interpret 6q 
as the exogenous part in the dynamics of pref and 9i as the intensity of the endogenous effects 
driving p^ef towards the current mid price level. 


3 Ergodicity 

In this section, we discuss ergodicity properties in our framework. To do so, we make some 
additional assumptions on the functions /*, Pi, u and d. 


Let 

Pt{x,A) := P[y(t) e ki|y(o) = x] 

be the transition probability at time t of a continuous or discrete-time Markov process Y with 
state space y. In this work, we say that the process Y is V-uniformly ergodic if there exists a 
coerciv^ function 1/ > 1, an invariant distribution vr, r G (0,1) and R > oo such that for any 
X E y and t G (or N’*' in discrete-time). 


\Ptix,.)-7r{.)\\v<Rr^V{x), 


(5) 


where we write ||.||v for the E-norm of a signed measure, see Meyn and Tweedie (1993, 2009). 


In continuous time, the main idea to prove such property is to design an appropriate Lyapunov 
function E: E —)■ (l,cx)), on which the following negative drift condition is satished for some 
7 > 0 and B > 0: 


SVfe) := Y.^vy'lVW-nv)] 

y'¥=y 

< --fV{y) + B. 


Then by Theorem 6.1 in Meyn and Tweedie (1993), the Markov process Y is non-explosive and 


V-uniformly ergodic. Furthermore, by Theorem 4.2 in Meyn and Tweedie (1993) it is positive 
Harris recurrent. Note that the same kind of method is used in Abergel and Jedidi (2011) in 
order to show ergodicity properties of Zero-intelligence models. 


As mentioned in the introduction, the LOB’s ergodicity implies here the existence of a unique 
invariant distribution for the state vector q. This is relevant for explaining the stylized empirical 
distribution of the LOB state. Mostly, as we will see in Section]^ the ergodicity analysis is the 
basis for proving the diffusive limit of the reference price process. 


3.1 When pref stays constant 

We hrst discuss the V-uniform ergodicity of the process q{t) when assuming u{q) = d{q) = 0 
in Equation (|^. Recall that the unused values of fi{q,n) and gi{q,n) in the dehnition of the 
queue dynamics, that is when q ± ncj ^ O, are set to zero. With the convention 0/0 = 0, we 


®|E(a;)| —)■ - 1-00 as ||a:|| —>■ oo. 
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define 


/;(?) 

■■= '^fi{g,n) 

n 


■= ^gi{g,n) 

k{q,n) 

fi{g,n) 

/;(?) 

ki{q,n) 

gi{g,n) 

9K9) 

Gk^^\z) 

:= ^z'^li{q,n) 

n=l 


00 

:= y^^z'^ki{q,n 

n=l 


Thus, when f*{q) > 0 (resp. g*{q) > 0), li{q,.) (resp. ki{q,.)) is a probability measure on 
N"*" with moment-generating function (resp. We make the four following 

assumptions. 

Assumption 3.1. For any order book state q and any i > ibestask, gi{q,n) = 0 for any n > qi 
and for any order book state q and any i < ibestbid, fi{q, n) = 0 for any n > —qi. 

Assumption 3.2. There exists z* > 1 such that for any q and i, < oo and G^’^'^{z*) < 

oo. Furthermore, there exists L > 0 such that for any i, 

z^l q 


Assumption 3.3. There exist r > 0 and U > 1 such that 

1 _ GSA<7(;^-1) 


lim sup [ff{q)-g*i{(l) r^fiqf^_ / J < 

_ X _ 

lim sup [g-iQ) - fiiQ) r 9 ,G(G-l ^ ^ 

{q,i):qi<-U,i<ibestbid ^ G) ^ 


( 6 ) 


Assumption 3.4. For any z > 1, 

B,(z) := 


inf {Gl-‘-‘‘(z) - 1) > 0 

(f?ib -Qi^^d^^bestask 

Bg{z) := inf {G^’^’^z) - 1) > 0. 

{q,i) '■qi'Zi U,i'^ibestbid 


To understand the practical meaning of these assumptions, let us consider the following example 
where the pure order book state jumps are assumed to have constant size equal to one. In 
such situation, the four assumptions above can be rewritten as follows and are much easier to 
interpret. 
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Example 2. LOB model with constant order size. 

• For n >2, fi{q, n) = gi{q, n) = 0 for any q E fl. 

• There exists L > 0 such that for any i G {—K, and q E Q, 

fi^Qj ^)'^i>ibsstbid diidll ^)'^i<ibestask ^ 


• There exist r > 0 and U > 1 such that 

sup [fi{q, 1) - gi{q, 1)] < -r 

sup [gi{q, 1) - fi{q, 1)] < -r. 


Basically, Assumption 3.1 says that a bid/ask limit cannot become an ask/bid limit in a single 
queue update event, that is the queue size cannot revert its sign in a single jumpj^ From 
Example we see that Assumption 3^ essentially states that the total intensity of the order 
insertion processes in the bid and ask side remains uniformly bounded with respect to the state 
of the LOB. Assumption 3.3, which is the most important for the system’s ergodicity, forces 
the individual queue sizes \qi\ to decrease when they become larger than a certain threshold. 
From these assumptions, we obtain the following theorem proved in appendix for the Markov 
process q{t). 


Theorem 3.1. When u = d = under Assumptions \2.1\ \3.1\ \3.2\\3.t^ and \3.4\ the continuous¬ 
time Markov jump process q{t) is non-explosive, V-uniformly ergodic and positive Harris recur¬ 
rent. 


Consider now the embedded Markov chain q{n) defined by q{n) = q{Jn), where Jn is the time 
of the n-th jump, and q{Jn) the state of the LOB after this event. The study of the embedded 
Markov chain is an important step in order to obtain the diffusivity of the price process in our 
setting. We write 


<(?) = 


f:iq) 


E^[ft{q) + 9n(l)] 


, bKq) = 


9iiQ) 


Ei[fn9)+9l{q)y 


for the proportions of queue size increases and decreases, and replace Assumption 3.3 by the 
following one. 


Assumption 3.5. There exist r > 0 and U > 1 such that 

lim sup [a*(g) - b*{q) . - —] 

_ l — nf,rq(z-0 

lim sup [b* (g) - a* (g) -—] 


< —r 

< —r. 


The following theorem is proved in appendix. 

Theorem 3.2. When u = d = 0, under Assumptions \2.1[ \3.1\ \3.^ \3.Ji\ and \3.^ the embedded 
discrete-time Markov chain q{n) is V-uniformly ergodic and positive Harris recurrent. 


^This assumption is not really mandatory but is realistic and technically quite convenient. 
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3.2 General case 


We are now interested in the case where u and d are no longer hxed to 0. Recall that q{n) 
represents the state of the LOB after the n-th event and Prefin) is the reference price (the 
center of the LOB) after the n-th event. We thus consider here the process of reference price 
increments c(n) (since the reference price itself is of course not ergodic), defined as the reference 
price change at the n-th event: 


c(n) = Pref{n) - Prefin - 1), 

and the embedded chain Yin) = (q'(n), c(n)), n G N, with c(0) an artihcial starting value. The 
process Yin) remains obviously Markovian. For some z > 1, let 

K 

r'([9.c])= 

i=—K,i^0 

We make two additional assumptions for the general case. 

Assumption 3.6. There exist z > 1 and L'^ > 0 such that for Qk, Q-k four 

random variables such that ~ ~ Qk ~ ttr and Q-k ~ '^-k- 

E[W([Q*’^", c])] + E[R^([Q'^“, c])] + E[ 2 l®^l-^] + E[zl'3-^l-^] < L^. 

Assumption 3.7. There exists a finite set W <zVL such that the upper bound of the proportion 
of reference price jumps in any order book state q is smaller than one on VL/W : 

_ ujq) + djq) _ 

geo/w l^iUiio) + Qiio)] + uiq) + d(g) 


dec 


Assumption 3T is technical and imposes some regularities on the four distributions tt*"", tt 
hk and ti-k- Assumption |3.7| ensures that a reference price change is not the only possible 


event except for a hnite number of LOB states. Under these assumptions, we have the following 
theorem proved in appendix on the ergodicity of the embedded chain Yin). 


Theorem 3.3. Under Assumptions \2.1\ \3.1[ \3.2\ \3.J\ |R5[ \3.(^ and |R?t the embedded discrete¬ 
time Markov chain U(n) = (g(n),c(n)) is V-uniformly ergodic and positive Harris recurrent. 


4 Scaling limits 

We are now interested in the scaling limit of the reference price process. Let Jj be the time of 
the Uth jump of the process. Let 

Nit) = inf{n, < t} 

be the number of events until time t, with the convention inf{0} = 0. Let 2’(n) be the 
cumulative price change until the n-th event, that is 2’(0) = 0 and for n > 1: 

n 

i=l 
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We have 


Z{N{t)) =pref{t) -pref{0). 

Thus it represents the reference price at time t recentered its starting value. We show in this 
section the diffusive behavior of Z{N{t)) as n tends to inhnity. 


Consider again the embedded chain Y{n) = {q{n),c{n)). From Theorem 3.3, Y{n) is V- 
uniformly ergodic towards an invariant distribution tt*. We have the following theorem for 
the rescaled price process in event time := 

Theorem 4.1. Under Assumptions 2.1. 3.1\ 3.1^ 3.4, 3.5, 3.6 and 3.1, z/E^*[c(0)] = 0, then 
the series 


cr^ = [cq] + 2 ^ E^. [coCn], 


(7) 


n=l 


converges absolutely. Furthermore, ifY{0) ~ n*, we have the following convergence in law in 
D[0, oo).- 

S^^\t) aB{t), 

where B{t) is a standard Brownian motion. 


Proof. This theorem is a direct application of Theorem 19.1 in 

Billingsley 

(2009). Indeed, the 

sequence c„ is clearly stationary and ergodic in the sense of 

Billingsley 

2009 

) (for example 


since it is stationary and mixing). Moreover, it has a hnite second order moment and for all n, 
E^. [Cn] = E^. [co] =0. □ 

Theorem 4.1| shows that in event time, the large scale limit of the reference price process is a 
Brownian motion. However, the most relevant question is that of the large scale limit of the 
reference price in calendar time. Thus we now consider the process 


\/n 


To prove the diffusivity of S^'^^t), we need a last assumption which is a bound of the expected 
value on the waiting time between two events. 

Assumption 4.1. There exists some m > 0, such that 

inf { J2^fiiq) + 9*iiq)) +u{q) + d{q)} > m. 


Let Tn be the inter-arrival time between the n-th and the (n — l)-th jumps of the Markov process 
X. We then have the following theorem proved in appendix for the long term behavior of the 
reference price in calendar time. 


Theorem 4.2. Under Assumptions 2^, 3^, 3^, 3^, 3_F. 3.(^ \3. 7| and 4jJ_> the process 


{q{n),c{n),T{n)) is positive Harris recurrent. Furthermore, z/E.^-* [c(0)] = 0 and T(0) 
then 

^ ^ B{t), 


TT 


V'E^*.[r(l)] 


with TT** is the invariant distribution of {q{n), c{n),T{n)) and a defined in ([^. 

Theorem |4^ discusses the scaling limit of the underlying reference price. However, the difference 
between this price and the more usual Pbestudit), Pbestaskif) or Pmidit) being bounded by 2K, 
the same result applies replacing the reference price by any of those prices. 
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5 Some specific models 

The Markovian setting proposed in this work allows us for a wide range of possibilities for 
modeling order book dynamics. The goal of this section is to give some natural and tractable 
examples of models, essentially already introduced in the literature, which can be seen as 
particular cases of our general framework. Together with the dynamics of the models, we 
provide sufficient conditions so that the assumptions made in the previous sections are satished 
in these specihc models. Thus the ergodicity and diffusive scaling properties apply in all these 
models. 


5.1 Best bid/best ask Poisson model (Cont and De Larrard (2013)) 


The basic idea of this hrst model, inspired by that introduced in Cont and De Larrard (2013), 
is to use a constant spread size (hxed to 1 tick) and to consider only two limits in the order 
book. 

Example 3. Poisson model with K = 1. 

• We take K = 1, 0^“™* = 1 and assume that the functions fi, gi, u and d have the following 
forms, with 0 < A < p < oo.- 


and satisfy Assumption\3.(^ and 


/i(g,u) 




/-i(g,u) 

/ilg_l<oln=l 

Wi(g,u) 


u{q) 

= 91q-^^=o 

d{q) 

= ^lg_i=0- 


for any q_i > 0, qi E h, 
for any qi < 0, q_i G Z, 

7r*-(g_i,gi)=7r''-(g_i,gi)=0. 

Note that here, the boundary distributions ttk and 7t_k are no longer needed, since the order 
book reinitialization probability is set to one. 


In this model, the role of pref is very close to that of Pmid, which splits the order book into 
two parts: the bid side (Q-i) and the ask side (Qi). The limit order insertion, cancellation 
and market order insertion processes are assumed to be independent Poisson processes. The 
size of these orders is assumed to be constant and Pref jumps with rate 0 when none of the 
queues Q±i is empty, with rate 9 to the right side when Qi is empty, with rate 9 to the left side 
when Q_i is empty. When the value of 9 is very large, the price jump is almost instantaneous 
as soon as one of the two queues becomes empty. In that case, this model becomes very close 
to that proposed in Cont and De Larrard (2013), where an inhnite rate is used (note that the 
convergence of the rescaled price process can still be proved with some minor modifications in 
such case of inhnite jump rate). 
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5.2 Poisson model with K > 1 


It is natural to try to extend the previous Poisson model in order to include more queues in 
the order book and to allow for a spread size different from one tick. In such model, the role of 
Pref is slightly different since it is not necessarily the mid price. Again, p^ef can be understood 
here as the underlying efficient price that determines the order arrival intensities at different 
price levels. Now buy/sell limit orders can be inserted both on the right side and on the left 
side of Pref- 

Example 4. Poisson model with K > 1. 

• The functions fi, Pi, u and d have the following forms, for i = —K ,..., K: 





f^'i'^i'^ibestbidi.Q') 




^ ( q ) 

fiiQ^n) 

= 0,/or n > 1 


= 0,/or n > 1 

u{q) 

^'^bestask 

d{q) 

_ n 

'^bestbid {q) * 


71^, n ^and 7r'^“ satisfy Assumption 3.6. 


• For any i,j G {—K, i < j, we have 


P—i ^ ^ 0 

0 < < Oj. 


Limit orders, cancellations and market orders (which consume the quantities at the best offer 
limits) are modeled by independent Poisson processes with different intensities ^ ^buy/seii 

and depending on the distance from their target price to Pref- We assume bid-ask 

symmetry in this model, that is = Ai®“, and thus we omit the 

index buy/sell in the above equations. Remark that the intensity of the buy/sell market order 
flow at the best limit y* is a function of i, that is the position of the best limit with respect 
to the reference price. This allows us to model the fact that market participants have different 
behaviors towards the best limit, depending on their evaluation of the reference price. 


The reference price jump dynamics is modeled by a sequence of increasing rates 9i, i G {—K, 

This means that the larger the index of the best ask queue, the larger the probability of Pref 
to increase and the smaller the index of the best bid queue, the larger the probability of Pref 
to decrease. Note that in this model, we no longer assume any specific value for the reinitial¬ 
ization probability 0^“™* and use Assumption 3T to impose some properties on the boundary 
distributions vr^ and 7t~^ and the initialization distributions and 7r‘^“. 


15 




5.3 Zero-intelligence model 


We now present a different way of extending the Poisson model with iP = 1 in order to include 
more queues in the order book. This modeling approach where two reference prices are used is 


called Zero-intelligence model and is introduced in Smith, Farmer, Gillemot, and Krishnamurthy 


(2003). It is also considered in Abergel and Jedidi (2011) and is the basis of Cont, Stoikov, 


and Talreja (2010). We dehne (j){i,j) as the absolute distance (in number of ticks) between the 
queue Qi and Qj ( 0 (i, j) = \i — j\) and make the three following assumptions. 

Example 5. Zero-intelligence model. 

• The functions fi, Qi, u and d have the following forms, for i = —K ,..., K: 


fiiQ, 1 ) = 

1 ) = 
fi{(l,n) = 

= 

u{q) = 
d{q) = 

7T^, 7T~^, and 7r'^“ satisfy Assumption 3.6 


^(p6,ibestbidlQ))'^'^^^bestbidiQ) — ibestbidlA \9.i\h‘<l>{i,'ibestaskiSl)')'^'^'^Uestbid{,Q) 

'^0(^!^6estasfc (q)) (?) (?) \9.'i.\T(l>{i^^bestbid{,A^'^'^'^^bestask{u) 

0,/or n > 1 
0,/or n > 1 

^^bestask (?) 

^ ^bestbid (?) * 


• For any i,j G {-K, ...,K], > 0 and > 0. 

• For any i,j G {—K, ...,K}, i < j, we have 

0<ei< Qj. 

In this model, Pref is no longer an underlying efficient price determining the order arrival 
intensities. These intensities now depend on the positions of two different prices: Pbestud and 
Pbestask- Limit Orders, cancellations and market orders are still described by independent Poisson 
processes. Buy/sell limit orders are inserted in the queues to the left/right side of the best 
ask/best bid price, with intensities depending on the distance between their price level and the 
best ask/best bid price {.^(f>{i,if,,,„t,,„k/bestbid))'^ cancellations of buy/sell orders are sent to the queues 
on the left/right side of the best ask/best bid price, with intensities being linear functions of 
the queue sizes i\qi\F<t>{i,ib^stask/b^stbid))'^ market buy/sell orders are sent to the best ask/best bid 
queue, with intensity 7 . The reference price p^e/ now provides the center of the 2K dimensional 


moving frame representing the LOB’s state and the same modeling approach as in Section 5.2 
is used for its dynamics. 


5.4 Queue-reactive model (Huang, Lehalle, and Rosenbaum (2015)) 

In Huang, Lehalle, and Rosenbaum ( 2015[ ), the Queue-reactive model for order books is intro¬ 
duced. This model takes into account the influence of the order book’s state in determining 
the order arrival intensities (in a much more general way than considering only the position of 
the best bid and best ask queues). The Queue-reactive model assumes that no buy/sell limit 
order can be inserted on the right/left side of pref and uses the following assumption instead of 
Equation (|^ for the dynamics of the jumps of Pref- 
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Assumption 5.1. Whenever pmid increases (resp. decreases), pref increases (resp. decreases) 
by a with probability 9, provided that = 0 (resp. g_i = Q) at that moment. Therefore, changes 
in the value of pref are possibly triggered by one of the three following events: 

• The insertion of a buy (resp. sell) limit order within the bid-ask spread while Qi (resp. 
Q-i is empty). 

• The cancellation of the last limit order at one of the best offer gueues. 

• A market order which consumes the last limit order at one of the best offer gueues. 


With some minor modifications in the proof of Theorem 4^, one can prove that the scaling limit 
of Pref in the Queue-reactive model is a Brownian motion. As explained above, in this model, 
changes of Pref are triggered by events that modify the mid price. Here we propose a slightly 
modihed version of the Queue-reactive modej^ by considering the following four assumptions 
(note that the state space H is reduced in that case to := {q G H, gili<o < 0, gjli>o > 0}). 

Example 6. Queue-reactive type model. 

• The functions f, gi, u and d have the following forms, for i = —K ,..., K: 


fiiQ, 1) 

— •^|i|(<?i)li>0 + h|i|(~5'i)li<0 


T hljl (O'i) li>0 

fiiqw) 

= 0,/or n > 1 


= 0,/or n > 1 

u{q) 

= 

d{q) 



with A|i| and p\i\ non-negative functions defined on N, with /i|j|(0) = 0. 
• We have 


sup (Aj(gj)) < L < oo. 

• There exist r > 0 and U > 1 such that for any qi > U and any i G {1, ..., K}: 

Ai(gi) - Piiqi) < -r. 


• There exists m' > 0 such that for any i G {!,..., iV}.- 

inf [Xfiqi) tifiqi)] > m'. 


Compared with Assumption 5.1, changes in pref are now driven by the relative position of the 
mid price in the current order book state. Nevertheless, we can see that the two approaches 
are quite similar. In this model, p^e/ always stays between Pbestask and Pbestud (since /i|i|(0) = 0 
implies that the queue sizes on the left/right side of p^e/ never become positive/negative). Such 


SModel I in|Huang, Lehalle, and Rosenbaum| (|2015[) is used to describe the queue dynamics during constant 


reference price periods. 
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model gives us a much larger choice on the intensity functions A|j| and ^\i\ than when assuming 
Poisson flows. Furthermore, with enough data points, these functions can be estimated in a 
non-parametric way, as done in Huang, Lehalle, and Rosenbaum (2015). Finally the state- 
dependent approach provides us very interesting insights about the way the order book state 
influences market participants decisions and the mechanism making the empirical distribution 
of the order book arise from these decisions, see Huang, Lehalle, and Rosenbaum (2015). 


6 Conclusion 


In this work, we extend the order book modeling approach proposed in Huang, Lehalle, and 


Rosenbaum (2015) to a more general Markovian framework, allowing to take into account most 


relevant features of LOB dynamics such as: 


• Dependencies between the order arrival processes and the LOB state. 

• Endogenous movements of the underlying efficient price and inffuence of the LOB state 
on its dynamics. 

• Exogenous movements of the underlying efficient price. 

• Randomness in the size of the orders. 


The ergodicity of the LOB system and the diffusive limit of the rescaled price process are es¬ 
tablished under general assumptions. Finally, to illustrate the usefulness and the relevance of 
our approach, several examples of classical models which can be seen as particular cases of our 
general framework are presented. 


To get a fully satisfying model, a last step would probably be to allow for a non-Markovian 
component in the market order flow (since the Markov assumption is probably quite reasonable 
for the limit order and cancellation flows). This can for example be done using self-exciting 
processes, as in Abergel and Jedidi ( 2015[ ). However, except for very specihc cases (exponential 
Hawkes processes for example), adding such non-Markovian component would certainly require 
revising completely the mathematical approach to the model. 


A cko wle d ge me nt s 

We thank Frederic Abergel for helpful discussions. 
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7 Appendix 


7.1 Proof of Theorem 3.1 


Let us denote by Q the infinitesimal generator matrix of q{t) when u = d = 0. The infinite 
matrix Q has the following form: 

Qq,q+nei — 

Qq,q—nei 

q£n,q'^q 

Qq,q' = 0, otherwise. 


For some 1 < z < z* (recall that 2 ;* is defined in Assumption 3.2), let us consider the function 

K 

V{q) = 

i=—K,i^0 

Since g* > 0 for i > ibestask, g* < 0 for i < ibestud and g* = 0 for i e {ibestbid, ibestask), we have 

QV{q) := 

Q'¥=g 

00 

= E - z-^^-^)U{q, n)] 

'^^'^bestbid ^=1 


E EK- 

'^^'^bestask ^—1 


qi+n-U _ ^qi-U 




+ E - z-^)f,{q,n) + (z"-^ - z-^)g.{q,n)]. 

'^^{'^bestbid^'^bestask') ^ ^ 


Then by Assumption 3.1 


av(.q) = 5 ^ -z-^--")gi(q,n) + (z 

'^^‘^bestbid ^^ 


-qi-n-U _ -qi-U 


Z * )fi(.q,n)] 




'^^'^bestask ^^ 


'^^{'^bestbidAbestask') ^ ^ 


. 5 : ^ 

'^^'^bestbid 


,n 


-*-‘'E(‘”-1)I9‘(«' 

n=l 
00 


i>i 


— ^bestask 


n=l 


Zn ^ 

9i{q,n) 

Z-n 


+ Y z + 9i{9,n)]. 

'^^{'^bestbidAbestask') ^ ^ 
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Using the definition of fi{q) and g*{q), we get 


QViq) = Y1 ^■*"‘'[9.*fe)(G»’‘’’(z)-l)-/;fe)(l-G'«(.--‘))] 

'^^'^bestbid 

'^^'^bestask 

+ E z-"[/-(9)(G«-'(2)-1)+9-(9)(G»«(.-)-1)]. 

(^hestbid ?^besfasfc ) 


Moreover, since for z > 1 we have < 1 and G^l^'^'\z) > 1, we obtain 

Qy{q) < 


E - 1)(9*(9) - fM) / ) 


G9-‘.«(;) - 1 


i-i'^ibestbidiQ.i'^ C/ 

+ E - 1)(/-(9) - 

i-i^ibestask 

+ E 9.*(9)G‘’«(9) + E A'(9)G'-‘-’(9) 

^•‘^^^bestbzd ^ '^•'^'^'^bestask ^ ^ 

+ E [/;(9)G"'»(9) + 9-(9)G"'«(9)| 

^G (ibestbid t'^bestask') 


< 


i'i^ibestbid^Qi^ C/ 




E 9''--‘'(G/«(9) - l)(/-(9) - 9;(9)^gj^^^) 


i’i^ibestask 

■ EW(9)G'’‘''(9)l.>i...,„ + 9.-(9)G‘''‘-’(.-)1.<... 


Now note that 




is an increasing function of 2 ;. Thus by Assumption 3^, we can find z' > \ such that for any 
z < z', q & VL and i G [—K, ..., A']: 

(/;(9)G-''‘’’(9)l,>i.„,„ + 9-(9)G9’«(9)li<,....„J < L. 

From Assumption |3.3[ we see that we can find some 5 with 1 < z < z' such that for any 
(g, i),qi< -U, i < ibestbid, for any 1 < z < z: 


9 . ( 9 ) - fi ( 9 ) a,_ 1 < -G 

and for any (g, i),qi> U,i> ibestask, for any 1 < z < z: 

fi ( 9 ) - 9 . ( 9 ) ^ < - 9 . 


( 8 ) 


(9) 
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Thus taking ^ in the dehnition of the function V satisfying 2Kz ^ > 1 (so that the function 
V is coercive) and 1 < z < z, we obtain 

QV(g) < -r z-‘^^-^(G^’^’'^(z)-l) 

-r - 1) + 2KL. 

i'-i^ibestask 

< -rBg{z) Y1 

'^•'^'^'^be.sthidiQi ^ 

-rBf{z) Y2 + 2KL. 

'i-'i^ibestask 

B := nim.{Bg{z), B f{z)) > 0. Therefore we get 


By Assumption 3.4 


QV{q) < -rB ^ z^^‘^-^ + 2KL 

r.\qi\>U 

< -rB^^I'?*l-^ + 2A:(L + r5) 

i 

= -rBV{q) + 2K{L + rB). 

Finally, remark that in our setting, any compact set included in hi is hnite. A singleton being a 


petite set and a hnite union of singletons remaining a petite set, see Proposition 5.5.5 in Meyn 


and Tweedie ( 

2009 

), we get that all the compact sets are petite. Therefore by Theorem 6.1 

in 

Meyn and Tweedie 

(1993), q{t) is non-explosive and V-uniformly ergodic. Furthermore, by 

Theorem 4.2 in Meyn and Tweedie (1993) it is positive Harris recurrent. 


7.2 


Proof of Theorem 


3.2 


For some 1 < 2 ; < 2 ;*, set again 


V{q)^ Y, 

i=—K,i^0 


We write Vq^q^ the transition probability from q to q'. In the same way as in the preceding 
proof, we have 


AV(q) := ^,.,.(^( 9 ') - l"(?)) 

g'eo 

'^^'^bestbid 

'^^'^best ask 

’^^i^^bestbidj'^bestask') 


Following the same method as in the proof of Theorem 3T, we can easily hnd 1 < z' < z* and 
B > 0 such that taking 2 ; = ^' in the dehnition of V, it is coercive and we get 


AV{q) < -rV{q) + B. 
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Now define the set C : = 
Furthermore, we have 


{q, rV{q) < 2B}, C is obviously a hnite set and is therefore petite. 

AV{q)<-'^V{q) + Bl,^c. 


Thus by Theorem 16.1.2 in Meyn and Tweedie (2009), q{n) is V-uniformly ergodic. 


Eventually, the fact that the chain is positive Harris recurrent is deduced from Theorem 9.1.8 
together with Theorem 15.0.1 in Meyn and Tweedie (2009). 


7.3 


Proof of Theorem 


3.3 


For ease of notation, we write V instead of V^. Let (5“^“, Qk, Q-k be four random 

variables such that ~ ~ Qk ~ t^k and Q-k ~ t^-k- We dehne 


u*{q) 

d*{q) 


n*{q) 

Er 

E_k 

IE 

lE^dec 


_M(g)_ 

Ei[/r(9) + + u{q) + d{q) 

_%)_ 

+ 9* {a)] + u{q) + d{q) 

EJ/;(g)+17?(g)] 

Ei [/*(?) + 9*i{.q)] + u{q) + d{q) 
E[z^^-^] 

E[H([Q*"^c])] 

E[H([Q"“,c])]. 


Remarking that R([g,c]) does not depend on c, we write from now on V{q) instead of R([g, c]). 
Moreover, we set V[q^c]\q',c'] as the transition probability from state [g, c] to state [g', c'] and Vq^q/ 
as the transition matrix of the embedded chain g(n) when u = d = 0. Using the form of the 
inhnitesimal generator Q, we deduce 


Ar(|,,c]) := 5^ PMj,vi(r(«') - r(«)) 

(g^c^)G^^X { —q,q:} 

= n-iq)Y,V,,^(V(q')-V{q)) 

q' 

+u*{q) [(1 - - U(g))] 

+d*(g) [(1 - - U(g))]. 


By Assumption 3.6 


we have 


AU([g,c]) < n*{q)Y,KAy{q)-y{q)) + ‘^L^+ ‘^L\ 

g'eo 


Then as in the proof of Theorem 3^, we can End 1 < 
z = z' in the dehnition of V, it is coercive and we get 


z' < z* and R > 0 such that taking 


AU(g) < -n*(g)rU(g)+ Rn*(g)+ 2L^ + 2L’^. 
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Moreover, Assumption |3.7 ensures that for any q except those belonging to the finite set W, 
n*{q) > M with M G (0,1]. Consequently, 

AC(g) < —Mrl/(g) + i? + + 2//^^ + z/l[g^c]eiVx{-o,a}; 

with 1 / > 0. Now define the set C := {[g, c], MrV (q) < 2B + AL^ + AL^}. Being finite, CVJW 
is a petite set and we have 


Mr 


Al/([g, c]) <- c]) + {B + 2L^ + 2U' + z^)l[q,c]ecuw- 


Hence from Theorem 16.1.2 in Meyn and Tweedie (2009), Y{n) is V-uniformly ergodic. 


Eventually, the fact that the chain is positive Harris recurrent is deduced from Theorem 9.1.8 


together with Theorem 15.0.1 in Meyn and Tweedie (2009). 


7.4 Proof of Theorem 4.2 


7.4.1 Preliminary lemma 

We start with the following preliminary lemma. 

Lemma 7.1. For the Markov chain {q{n),c{n),T{n)), the Cartesian product of any finite set 
included in fl x {—a, a} and is petite. 

Proof. We first show that for any q E fl and c G {—a, a}, the set g x c x IR+ is petite (actually 
small). We define the measure r'g^c,T which is so that for any B G B{fl x {—a, a} x M^), Ug^criB) 
is the transition probability from [g, c, r] to H in a single step: 

^q,c,T(7?) 

Recall that in our framework, the transition probabilities from {q{n),c{n),T{n)) depend only 
on the value of g(n). So we can write Ug^criB) as I'qi^B). In the sense of Equation (5.43) in 
Meyn and Tweedie (2009), the transition probability V[g^r,T],B can be seen as a sampling kernel 
for the Markov chain (g(n), c(n), r(u)), using the Dirac measure at point 1 on IC as sampling 
measure. Moreover, for any r G and any B G B{fl x {—a, a} x M"*"), we have 

Big, c,t],b > J^q{B). 

Therefore the set q x cx M’*' is petite. Then, as the union of two petite sets remains petite, see 
Proposition 5.5.5 in Meyn and Tweedie (2009), we have the result. □ 


7.4.2 A law of large numbers 

In the next proposition, we give a law of large numbers for the inter-arrival times, which is the 
key element to establish the diffusive behavior of the price in calendar time. Within the proof 
of this proposition, we show that {q{n),c{n),T{n)) is positive Harris recurrent. 

Proposition 7.1. Let Ti be the inter-arrival time between the i-th and the i — 1-th jumps of the 
Markov process X. Under Assumptions \3.1\ \3.1\ \3.f \ |R5[ \3.(^ \3. 7| and 4jJ_, almost surely, we 
have 

1 

^ E^*.[r(l)], 

* ^ n.—^ro 


n 


2 = 1 


with 71** the invariant distribution of {q{n),c{n),T{n)). 
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Proof. First we show that the Markov chain {q{n),c{n),T{n)) is positive Harris recurrent. For 
the Markov chain {q{n), c{n)), we have already proved that a coercive function V can be found, 
such that the following drift condition is satisfied for some a > 0 and L < oo: 


AV{[q,c\)<-aV{[q,c]) + L. 

Now, for the Markov chain {q{n),c{n),T{n)), take 

y*{[Q,c,t]) = V{[q,c])+t. 

With obvious notation, we have 

AW([g,c,f]) = AV{[q,c])+E[q^c,t][r]-t 

< —aV{[q,c\) + L + Kx[T]—t. 

Taking a' = min{a, 1) and L' = L + 1/m, using Assumption |4.1 , we get 

AV*{[q,c,t]) < -aV{[q,c\) + L+ l/m-t 

< -a[V{[q,c])+t] +L+ l/m 

< -aV*{[q,c,t]) + L'. 

Now let C = {{[q,c,f\), a'V*{[q,c,f\) < 2L'}. We have 


AV*{[q,c,t]) < --W([g,c,f])+ L'l[g,c,t]ec- 


According to Lemma 7.1, the set C is petite. Thus we can apply Theorem 15.0.1 in Meyn and 


Tweedie (2009) to deduce that the Markov chain {q{n),r{n),T{n)) is positive recurrent and 


thus admits an invariant measure. 


Now remark that the function V* is unbounded off petite sets (using Lemma 7.1 together with 


the fact that any subset of a petite set is itself petite). Consequently, Theorem 9.1.8 in Meyn 


and Tweedie (2009) enables us to obtain that the Markov chain (g(n), c(n), r(n)) is Harris re¬ 
current. Therefore it is positive Harris recurrent. 


We end the proof thanks to Theorem 17.0.1 from Meyn and Tweedie (2009). 


□ 


7.4.3 End of the proof of Theorem |4.2 

We have 


_ Z([i»/E,..[r(l)|J) ^ . Z(N(nt)) _ Z(Ltn/E,..|r(l)]J) 

\/n ^ \/n \/n 


According to Proposition 7.1, the sequence of processes N{nt)/n converges to [r(l)]. 

Moreover, the limit of Z([fn/E^**[r(l)]J)/-y/n is continuous. Thus, using Skorohod representa¬ 
tion theorem together with continuity properties in Skorohod topology, see Proposition VI.2.1 
in Jacod and Shiryaev (2003), we get that the second term on the right hand side of the above 


equality tends to zero. Finally, from Theorem |4.1[ we get 

Z([fn/E^..[r(l)]J) 


n—^oo 

-)■ 


a 


n 


\/E^*. [r(l)] 


Bt. 


Combining these two results, we obtain the weak convergence of the rescaled price process to a 
Brownian motion with variance (J^/Ett** [r(l)], which concludes the proof. 
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